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The cross section of the 16 O(«,7) 20 Ne capture reaction is analyzed at low energies where the 
direct capture mechanism is dominant. For temperatures below Tg = 0.2 the resulting astrophysical 
reaction rate is about a factor of two higher than in a recent compilation whereas the energy 
dependence of the astrophysical S-factor and the branching ratios to the 20 Ne bound states are very 
similar to previous calculations. The validity of the widely used detailed balance theorem for the 
inverse 20 Ne(7,a) 16 O photodisintegration rate is confirmed for the special case of high-lying first 
excited states. 
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I. INTRODUCTION 



The astrophysical reaction rate of the 16 O(a,7) 20 Ne reaction is dominated by a series of resonances in a broad 
range of high temperatures above 0.2 x 10 9 K (Tg > 0.2). However, at low temperatures (Tg < 0.2) the direct capture 
O^l mechanism becomes dominant p). 

The 16 O(a,7) 20 Ne reaction is an important reaction during helium burning in stars at typical temperatures around 
Tg = 0.2 which is the temperature range where both direct capture and resonances are important. The 16 O(a,7) 20 Ne 

■ reaction blocks the reaction chain 3a -* 12 C(a,7) 16 O(a,7) 20 Ne because there are no resonances in the so-called Gamow 

window which is located at E w 1.22 (Z 2 Zf A rod T 6 2 ) 1/3 keV with a width of A = 0.749 (Z 2 Z|A red T 6 5 ) 1/6 keV; Z x 
(Ai) and Zi (A%) are the charge (mass) numbers of projectile and target, and A roc i = AiA2/(A x + A2) is the reduced 
mass number. At Tg = 0.2 one finds E — 390 keV and A = 190 keV. In several previous papers the S-factor at 

■ Eq = 300 keV is reported. The 2~ state at E x — 4967 keV in 20 Ne does not appear as resonance at E = 237 keV in 
[ the 16 O(a,7) 20 Ne reaction because of its unnatural parity. 

The inverse 20 Ne(r,a) 16 O reaction plays an important role in neon burning at significantly higher temperatures 
around Tg = 1 — 2 ,2j. The reaction rate of photodisintegration reactions is usually derived from the capture rate 
using the detailed balance theorem (see e.g. p|). The validity of this theorem for the case of the 16 O(a,7) 20 Ne and 



20 Ne(7,a) 16 O reactions with high-lying first excited states is analyzed in the last section IVT1 



Direct capture in the 16 O(a,7) 20 Ne capture reaction has been analyzed theoretically by 0, 0, IS IE S • Exper- 
imentally a non-resonant S-factor of S — 0.26 ± 0.07 MeVb for the ground state transition has been reported at 
energies around E = 2MeV 0, and a total S-factor of 5(300 keV) < 0.7MeVb was deduced, but this result has been 
questioned by Q. Later, an upper limit of S < 3.37 MeVb was reported which was derived by several data points 
between 1.4 MeV and 2.4 MeV; no extrapolation to 5(300 keV) is given in jlOj because of the contradicting energy 
dependencies of the S-factor in [3,01 . The theoretical predictions vary between 0.7 MeVb and 2.5 MeVb at 300 keV 

This work gives a new prediction for the direct capture cross section of the 16 O(a,7) 20 Ne reaction at low energies. 
The calculation is based on a systematic folding potential which is able to reproduce simultaneously scattering data for 
16 0(a,a) 16 and bound state properties of the system 20 Ne = 16 ® a jlj. The direct capture model is presented in 
Sect.[n] Relevant bound state properties are listed in Sect. IIIII The direct capture cross section is shown in Sect. II VI 
and the results are discussed in Sect.[V] After a brief analysis of the inverse 20 Ne(7,a) 16 O photodisintegration reaction 
in Sect. IV1I conclusions are drawn in Sect. IV 111 
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II. DIRECT CAPTURE, B{E/MC) VALUES, AND FOLDING POTENTIALS 



A. Direct Capture 

The direct capture (DC) cross section for El, E2, and Ml transitions is given in detail in Refs. [l2lll^|. Here only 
the essential properties of the model are repeated. The DC cross section is proportional to the square of the overlap 
integral of the scattering wave function xi > the electromagnetic transition operator Q E / MC ) and the bound state wave 
function uatl : 



<Jdc{E) 
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u NL {r)0 E ' MC X i(r,E) 
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The radial parts of the dominating El, E2, and Ml electromagnetic transition operators are given by: 
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(2.2) 
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with the so-called long wavelength approximation for p = k^r <C 1. From Eq. 1)2. 4|) and from the orthogonality of the 
wave functions it is obvious that the Ml contributions will be very small in the DC model (sec Sect. IIV(I . The exact 
form of the multipole operators in Eqs. l|2.2|l . (|2.3|l . and l|2.4f) was used in the following calculations. 
The coefficient of proportionality in Eq. (|2.1|l for El transitions contains a factor 



6(251) = [C(E1)Y 



Ai 



Z2 
A 2 



(2.5) 



which describes the isospin suppression of El transitions in N = Z nuclei between T = states. If one uses integer 
mass numbers, the El cross section exactly vanishes; using precise masses, one calculates very small El cross sections. 
In reality, small T — 1 isospin admixtures in the scattering and bound state wave functions may lead to a significant 
enhancement of the capture cross section compared to the DC model calculation [j| . 



B. B(E/MC) values 



Whereas Eq. (|2.1|l determines the overlap integral at an arbitrary energy E of the scattering wave function, the 
following Eq. 1)2. 6f) gives a similar ove rlap integral between bound state wave functions at their defined energies E Xt i 
and E x j for transitions Li — > Lf [l4lll5| : 



B(EC, Li — » Lf) = 



e^(2L / + l)(2T + l) 



with 



47r(2Li + 1) 



(LfOCtyLA) 2 



UN f L f (r)r c u NiLi (r) dr 



A^Z X + {-l) c AfZ 2 
(Ax+A 2 )£ 



(2.6) 



(2.7) 



The above formalism is also valid for quasi-bound states with E > which appear as resonances in the capture 
reaction. 

Obviously, one finds the same suppression of El transitions between T — states in N = Z nuclei because the 
suppression factors in Eq. H2.5[) and (|2.7|l are identical: C(E1) = (3\. This identity gives the chance to renormalize 
the DC calculation by the ratio of experimental and calculated B(EC) values 



Bcxp(EC) 



Bcalc(EC) 

leading to a theoretically predicted capture cross section a^: 

o-th(E) = v be x ctdc(E) 



(2. 



(2.9) 
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One expects ratios tbe ~ 1 for E2 transitions and tbe 3> 1 for El transitions. Because Ml transitions play only 
a very minor role (see Sect. HV|) . t^m = 1 was used in all cases. The expectations for the tbe values are confirmed 
by the calculations in Sect. IIIII 



C. Folding Potentials 

The essential ingredients for the calculation of the capture cross section and the B(EC) values are the potentials 
which are required for the calculation of the wave functions in Eqs. i|2.1|l and (|2.6|) . Based on the double- folding 
approach, a potential has been derived in ^l| which is capable of describing simultaneously elastic 16 0(a,a) 16 
scattering at energies above and below the Coulomb barrier and bound state properties of the system 20 Ne = 16 O <g> 
a. 

The potential V(r) is given by 

V(r) = V c (r) + XV F (r) (2.10) 

with the Coulomb potential Vc (r) of a homogeneously charged sphere with radius Rc and the folding potential Vf (r) 
which is scaled by a strength parameter A of the order of A w 1.2 — 1.3. 

V F (r)=J J Pi(ri) P 2(r2)v cS (E)d 3 ri d 3 r 2 (2.11) 

pi,2 are the densities of the a particle and 16 which are derived from electron scattering data [Tfij . and v c g is an 
effective interaction ^2 E ■ The Coulomb radius Rc has been chosen identical to the root-mean-square radius of 
the folding potential (r rms = 3.603 fm). 



III. BOUND AND QUASI-BOUND STATE PROPERTIES OF 20 Ne = 1b O ® a 

The wave function uml{^) which describes the relative motion of the a- nucleus system is characterized by the node 
number TV and the orbital angular momentum L. These values are related to the corresponding quantum numbers rii 
and li of the four nucleons forming the a cluster in the sd-shell: 

4 4 

Q = 2N + L = J2( 2n i + k) = X)* t 3 - 1 ) 

i=i i=i 

One expects Q — 8 for the ground state band and Q = 9 for the first negative parity band. The following procedure 
has been used for the calculation of the wave functions and their properties. First, the strength parameter A of the 
folding potential was adjusted to reproduce the energy of the state under consideration. Second, for resonances (which 
are quasi-bound states with E > 0) the theoretical width r^ h was derived from scattering phase shifts 6l by 



-ith 



(3.2) 



a (dS L /dE)\ E=EiQBtSL ^ /2 

The results are listed in Table [I] 

One finds that the strength parameter A is constant within about ±2 % for the ground state band with positive 
parity and within ±5 % for the negative parity band (or ±2 % if one neglects the 9~ state with its tentative assignment 
to this band). The average values A CV cn for the positive parity band are slightly smaller than the average A dd for 
the negative parity band. Additionally, a remarkable agreement is found between the calculated widths r^ h and the 
experimental widths T° xp especially for the negative parity band; for the ground state band one finds an overestimation 
of the order of a factor of 6.6 (14.5) for the 6 + (8 + ) state. 

It is interesting to note that the potential for the L = partial wave roughly reproduces also the properties of 
the 0|" state in 20 Ne at E x w 8700 keV. Using exactly the same potential as for the ground state (A = 1.2413) one 
obtains a broad resonance with Q = 10 at i?* h w 8080 keV with a width of several MeV. This compares well with 
the experimental values of the excitation energy of E* xp ~ 8700 keV and of the experimental width T° xp > 800 keV. 
However, it is difficult to define the excitation energy of this broad resonance from the potential model calculation 
because the calculated phase shift does not reach So = tt/2. Instead of Eq. 13.2f> a revised definition for broad 
resonances was used here: 

ti 2 

a {d8 L /dE)\ 
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TABLE I: Properties of the Q — 8 and Q = 9 states in 20 Ne. Experimental values have been taken from [Tg|. There is no clear 
assignment of 7~ and 9~ states to the 0~ band in [Tflj : the lowest tentative assignment was used. Minor numerical differences 
to Ref. [Till are a consequence of a different Coulomb radius Rc- 
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"derived from F a = 310 ± 30keV and r ao /r a = 0.51 ± 0.03 
'derived from T a = 220 ± 25keV and r ao /T a = 0.24 ± 0.01 



TABLE II: B(EC) values for electromagnetic transitions in 20 Ne. Experimental values have been taken from [Tsj|. 
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"vbe assumed as for 1 — > + ; no significant relevance for the 16 O(a,7) 20 Ne cross section, see Fig. HI 

In the same way, the properties of the broad 2 + state in the 0^ band are reproduced: using A = 1.2239, one finds 
Ef = 9130 keV and F^ h w 5MeV which compares to the experimental values = 9000 keV and T° x p w 800 keV. 

Using the potential strength parameters A as given in Tabled one may calculate B(EC) values using Eq. (|2.6|l and 
the formalism in Sect. Ill 3\ without free parameters. The results are listed in Table ITTI which includes also the ratio 
tee between the experimental and theoretical values as defined in Eq. I|2.8|l . 

The above simple two-body model reproduces the partial widths T a of the states under consideration with reasonable 
accuracy (see Table HJ, and the calculated reduced transition strengths for E2 transitions are in good agreement with 
the experimental values (see Table IT!]) . A prerequisite for the successful application of a simple two-body model is 
that the relevant states in 20 Ne have a dominating 16 ® a cluster structure. Such a structure can be expected from 
the double shell closure at 16 O, and it is confirmed by large spectroscopic factors close to unity which have been 
extracted from the analysis of transfer reactions. However, it is difficult to specify precise values for the spectroscopic 
factors as can be seen from the results given in compilations (3 Eil HH El and from recent experiments (see e.g. 

mm mm Hi). 
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IV. RESULTS 



A. The astrophysical S-factor 



The calculation of the capture cross section of the 16 O(a,7) 20 Ne capture reaction is now straightforward. The 
folding potentials (see Sect. Ill 0|) have been adjusted in strength to reproduce the bound state energies (see Sect. lIIIfl . 
Consequently, the wave functions of the DC model (see Sect. Ill A|) are well-defined. The absolute value of the capture 
cross section is adjusted using the ratios between the experimental and calculated B(EC) values (see Eq. (|2.8[1 and 
Sect. lITB)) . The results are shown in Figs. Q El and 01 




0.0 0.5 1.0 

E (MeV) 



FIG. 1: S-factor of the DC cross section of the 16 O(a,7) 20 Ne capture reaction. The sum of all transitions is shown as full 
line. Transitions to the + ground state (2 + first excited state; 4 + second excited state) in 20 Ne are shown as dotted (dashed; 
dash-dotted) lines. The Ml transition from the incoming L = 4 partial wave to the = 4 + bound state has a S-factor of 
less than 6 x 10 _6 keVb (off scale). The S-factor of the narrow 3~ resonance at E — 891 keV has been calculated using the 
Breit-Wigner formula (dashed line); this resonance is not included in the model space, see Sect. IVI 



Fig- summarizes the calculations. All possible El, E2, and Ml transitions from the incoming partial waves with 
< L < 6 to the three bound states with J 71 " = + , 2 + , and 4 + have been taken into account. Below E — 500 keV, 
about two third of the total cross section comes from the incoming L = partial wave to the first excited 2 + state. 
Further significant contributions from the incoming L = 2 partial wave to the ground state and to the first excited 
state are of the order of about 20% and 10%. The DC cross section below E = 500 keV is almost entirely given by 
E2 transitions. 

At energies above E — 500 keV the tail of the 1~ resonance at E = 1058 keV becomes more and more important. 
Although the width of this resonance is slightly overestimated by the potential model calculation (see Table HJ , the 
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adjustment of the B{EC) values according to Eq. 1|2.8[1 assures that the calculated resonance strength is in agreement 
with the experimental value. 

Transitions from partial waves with angular momentum L > 2 are practically not relevant for the 16 O(cv,7) 20 Ne 
cross section below 1 MeV. These transitions are suppressed because of the centrifugal barrier and because of the 
reduced transition energy for all transitions to the second excited state with J 71 ' = 4 + at E x = 4248 keV. The tail of 
the 3~ resonance at E — 2426 keV is not visible at energies below E = 1 MeV. 

The branching ratios to the ground states with J 71 ' = + and the first (second) excited state with J n = 2+ (4+) 
at E x = 1634 keV (4248 keV) are shown in Fig. El The branching ratios are almost independent of energy below 
E = 800 keV with a strong 78% contribution to the 2+ state and a 22% contribution to the + ground state. The 
contribution of the second excited 4 + state remains below 1 % in the whole energy range under consideration. 
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0.0 0.5 1.0 
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FIG. 2: Branching ratios of the DC cross section of the 16 O(a,7) 20 Ne capture reaction leading to the + ground state (dotted 
line), 2 + state at E x — 1634 keV (dashed line), and 4 + state at E x = 4248 keV (dash-dotted line; below 1%, practically not 
visible) . 



Experimentally upper limits have been determined for the S-factor to the + ground state (S — > + < 920keVb) 
and to the 2 + state (S — > 2 + < 2450 keVb) 10] from measurements between resonances at higher energies. The 
predictions from the DC model calculation are S -> 0+ = 802keVb and S -» 2+ = 2877 keVb at E = 300 keV. 
Because of the energy dependence of the S-factor, the calculated values do not contradict the experimental limits of 
In a similar way as in HI a total S-factor of S"(300keV) = 700 ± 300keVb is derived in |9j. However, the result 
of H has been questioned by |4| : the experimental data of |j at energies around 2 MeV were extrapolated to 300 keV 
using the probably unrealistic energy dependence of the S-factor from (see discussion in 0,0); it is stated that 
this result "should at least be multiplied by a factor of two" Q . Additionally, the upper limit of Q was derived close 
to the + resonance at 1991 keV; consequently the derived upper limit depends sensitively on the chosen resonance 
parameters (strength and width) of this resonance H3 . 

An estimate of the theoretical uncertainties can be done in the following way. Obviously, the calculated cross 
sections, which have been adjusted to experimental B{EC) values according to Eq. (|2.8I) . depend on the experimental 
B(EC) values or T 7 partial widths. Because T a 3> I\ for all resonances under consideration, the radiation width 
T 7 may also be obtained from the resonance strength tuj of a resonance in the 16 O(a,7) 20 Ne capture reaction. 
Typical experimental uncertainties for the resonance strengths are of the order of about 15% pj, and results from 
different experiments also agree with each other at the 15% level. Combining these uncertainties, this leads to a total 
uncertainty of about 20 % for the experimental B(EC) values or radiation widths T 7 . A second source of uncertainties 
comes from the adjustment of the potential strength to the binding energies of (quasi-)bound states (see TableCJ. The 
numerical adjustment itself can be done with negligible uncertainties; however, the model assumes a pure two-body 
configuration for all states, and consequently, the adjustment of the potential may be somewhat uncertain. But it 
has ben shown in H3 > that the potential reproduces experimental phase shifts with reasonable quality. Additionally, 
almost the same potential strength parameter A has been found for all states with positive parity; an uncertainty of 
about 2 % can be seen from Tabled The dependence of the S-factor or capture cross section on the potential strength 
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was calculated using a variation of the parameter A of ±2 %. For the dominating transition from the incoming L = 
partial wave to the 2 + bound state one obtains a variation of the S-factor of +14% (—15%) for a 2% enhanced 
(reduced) potential strength. Combining all the above uncertainties, a careful estimate of the total uncertainty of the 
S-factor is about 30%. 

B. The astrophysical reaction rate Na (p v) 

For the calculation of the astrophysical reaction rate Na (crv) in the temperature range where DC is dominating, 
one needs the S-factor from very low energies to about 500 keV. A second-order polynom has been fitted to the 
calculated S-factor which is shown in Fig. in linear scale: 

S(E) = 3808 keVb- 864 keVb x E/MeV + 1527keVb x (£/MeV) 2 (4.1) 




0.0 0.5 1.0 

E (MeV) 

FIG. 3: S-factor of the 16 O(a,7) 20 Ne reaction in linear scale (same as Fig. 0. For extrapolation to E — > 0, the S-factor has 
been fitted using a second-order polynom S{E) = 3808 keV b- 864 keV b x E/MeV + 1527 keV b x (E/MeV) 2 up to E = 0.5 MeV 
(long-dashed line). The S-factor of the narrow 3~ resonance at E — 891 keV has been calculated using the Breit-Wigner formula 
(dashed line); this resonance is not included in the model space, see Sect.fvl 

Because of the weak energy dependence of the S-factor the reaction rate can be calculated using the approximation 

N - = na G) 1/2 s{Ea) i 1 + ifc) exp ( ~ r) (42) 

with Eq and A as defined in Sect. U and r = 3Eo/(kT). Compared to 0], a somewhat larger reaction rate is obtained 
because of the larger S-factor at the relevant energies (see Fig.0J. 

At temperatures below Tg « 0.2 the reaction rate is almost entirely given by the DC contribution. Because of the 
increased S-factor at low energies, the reaction rate is enhanced by the same factor. The almost constant ratio below 
Tg = 0.2 in Fig. 0] arises from the fact that the energy dependence of the S-factor is practically the same in this work 
and in Q which was adopted from (see also Sect.EJ. At higher energies the reaction rate is dominated by the 
resonant contributions, and hence the increased DC cross section has no significant influence on the reaction rate. 

No distinction is made here between the reaction rates Na («}* under stellar conditions and Na (°" l; }i a b m the 
laboratory. Both rates are practically identical up to temperatures of about Tg = 3 There is not significant 
thermal population of excited states in 16 O which are located at excitation energies above E x — 6 MeV. The stellar 
reaction rate Na can directly be derived from experimental data. 
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FIG. 4: Ratio between the reaction rate Na (<? v) calculated from the DC model and folding potentials (this work) and the 
NACRE compilation lj. Because of the increased DC cross section at low energies, the reaction rate is about a factor of two 
larger than the NACRE rate below T9 = 0.2 where DC is dominating. 



V. DISCUSSION 



Because of the properties of the doubly-magic 4 He and 16 O nuclei a simple two-body model should be able to 
describe scattering phase shifts, bound state properties of 20 Ne, and the capture cross section of the 16 O(a,7) 20 Ne 
reaction with reasonable accuracy 0, 0, . The combination of the model with systematic folding potentials ensures 
a good description of elastic scattering cross sections and phase shifts at low energies. 

The strength parameter A is adjusted to the binding energies of low-lying states in 20 Ne; this guarantees the correct 
asymptotic behavior of the bound state wave functions at large distances which is important for the calculation of the 
16 O(a,7) 20 Ne capture cross section. The correct behavior of the wave functions at small distances, i.e. the number of 
nodes of the wave function and their positions, is obtained from Eq. (|3.1|) which avoids Pauli-forbidden states. The 
precise influence of antisymmetrization of the 20-nucleon wave function has been discussed in detail in |6( ; it has been 
pointed out "the smallness of the residual antisymmetrization terms" because of the shell closures in 4 He and 16 0. 

The potential in the chosen model automatically contains eigenstates at energies below and above the 16 0-a 
threshold. The quasi-bound states above the threshold appear as resonances in the 16 0(a,a) 16 scattering phase 
shifts and in the 16 O(a,7) 20 Ne capture cross section. As already pointed out in |(J, the term "direct capture" may be 
misleading in this case. Further details on such so-called potential resonances are given in [27l I2H I29I ] . However, the 
model contains only a limited number of eigenstates. It is not possible to include all resonances of the 16 O(a,7) 20 Ne 
reaction in this simple two-body model. E.g., the states with dominating 16 0(3~) ® a structure at E x = 4967keV 
(E = 237 keV; no resonance because of the unnatural 2~ parity) and 5621 keV (E — 891 keV) are not included in the 
model; the inclusion of such states requires e.g. a multicluster model 3j. The influence of the 3~ resonance at 891 keV 
at low energies is negligible as can be seen from Figs. ^ and |3J The S-factor of this resonance has been calculated 
from the Breit-Wigner formula using the known widths and strength [lj, Il8| . 

The calculation of the overlap integrals in Eq. (|2.1|l depends sensitively on the shape of the wave functions especially 
for wave functions with several nodes. Minor shifts of the nodes may reduce or enhance the overlap integral significantly 
because of cancellation effects between positive and negative regions of the integrand in Eq. (|2.1|l (see Figs. 2-5 of 
Ref. 0). This leads to a relatively sensitive dependence of the 16 O(a,7) 20 Ne cross section on the chosen potential, i.e. 
the potential strength parameter A. A variation of A by its uncertainty of ±2 % (see TableQJ changes the capture cross 
section already by about 15 % which has to be considered as the intrinsic model uncertainty. Further uncertainties as 
discussed in Sect. IIV Al lead to a total uncertainty of the calculated S-factor of about 30%. 

The adjustment of the DC cross sections of the 16 O(a,7) 20 Ne reaction using the ratios between experimental and 
theoretical B(EC) according to Eq. (|2.8|) gives the chance to include El transitions in the DC model. As stated 
earlier (see Eqs. (|2.5|) . (|2.7|l . and Sect.[nJ|, El transitions are suppressed in the DC model. A detailed microscopic 
investigation of El transitions in the 16 O(a,7) 20 Ne reaction was able to reproduce the approximate strength of the 
transitions, but showed problems with the branching ratios [jj. Therefore the usage of Eq. I|2.8|l seems to be the best 
way to predict the El part of the 16 O(a,7) 20 Ne cross section at low energies: the energy dependence of the cross 
section is provided by the model, and the absolute strength and the branching ratio are obtained from the adjustment 
to experimental data. The total contribution of El transitions to the 16 O(a,7) 20 Ne cross section remains small as 
long as energies significantly below the 1~ resonance at 1058 keV are considered. 

The final result for the total S-factor at 300 keV is 5(300 keV) = 3.7 ± 1.2MeVb which is slightly higher than 
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previous calculations 0,SSJESI3 but overlaps with the range of the results of Jj| 0, 0, @] within the estimated 
uncertainties. The result of [8| is significantly lower but has been questioned in The calculated slope a of the 
S-factor at low energies, usually defined as 



f fdS 



S(0) \dE 



(5.1) 



is in good agreement with previous calculations: a = — 0.23 MeV^ 1 from this work for the total S-factor is almost 
identical to a = — 0.21 MeV -1 for the dominating L = — * 2 + transition in |6j. The S-factor decreases slightly with 
increasing energy in all calculations 0, El El El ■ The only exception is the result of who finds an increasing 
S-factor with increasing energy with values of a « +0.6 MeV -1 . A possible shortcoming of Q is discussed in 
Branching ratios b, here often defined as 

b=^±^ (5.2) 
jo 

can be estimated from Fig. 5 of Q to be of the order of 3 < b < 3.5, and b rs 4.3 is given in In this work b = 4.5 
is found. Again this calculation agrees nicely with the microscopic calcaulations by |3J, |4jj. 

Any comparison with experimental data is difficult. No experimental data are available at energies below lMeV, 
and the extraction of the DC contribution from experimental data at higher energies is hampered by the dominating 
resonant contributions. As pointed out in the extraction of the DC cross section depends sensitively on the chosen 
resonance parameters, because the DC cross section is usually determined as the difference between the measured 
cross section and the calculated resonant cross section in the tail of a resonance with Breit-Wigner shape. As pointed 
out already in Sect. II V ~K\ the new calculation does not contradict the upper limits of [Tof . 

Finally, the astrophysical reaction rate Na (cr v) is about a factor of two higher than in the NACRE compilation I] 
below Tg = 0.2 because the S-factor at low energies is about 3.7 MeV b in this work compared to the adopted value of 
2 MeV b in [lj . A detailed analysis of the astrophysical consequences is beyond the scope of this paper; in general, the 
16 O(a,7) 20 Ne reaction rate is still low enough to block the reaction chain 3a ^ 12 C(a,7) 16 O(a,7) 20 Ne. Using the 
NACRE rate for the 12 C(a,7) 16 reaction and the new rate for the 16 O(a,7) 20 Ne reaction, one finds at Tg = 0.2 a 
ratio of about 5 x 10~ 5 between the 16 O(a 1 7) 20 Ne and the 12 C(a,7) 16 rates. For comparison, abundances calculated 
from the NACRE rates can be found in 



0£/) 



VI. THE INVERSE 20 Ne( 7 ,a) 16 O REACTION 

Reverse reaction rates A* under stellar conditions are usually derived from the detailed balance theorem 0, l3l| 
which reads for the case of the 16 O(a, 7 ) 20 Ne and 20 Ne(7,a) 16 O reactions 

exp -— (6.1) 



(av) \2irh J J V kT - 

because the partition functions do not deviate significantly from unity for the nuclei 4 He, 16 0, and 20 Ne up to 
temperatures of Tg ss 3 [1| because of the large excitation energies of the first excited states ( 16 0: E x = 6049 keV, 
0+; 20 Ne: E x = 1634 keV, 2+). The reaction Q- value is Q = +4730 keV. It has been shown in [32| that the detailed 
balance theorem remains valid for resonant transitions although the thermal occupation probability is extremely small 
for the first excited state in 20 Ne and negligible for 16 for all astrophysically relevant temperatures. 

The same arguments as given in |32) hold for the non-resonant case at temperatures below Tg = 0.2 (although the 
photodisintegration rate is practically negligible at such low temperatures). The reaction rate of the 20 Ne(7,a) 16 O 
photodisintegration reaction for a given state i in 20 Ne is given by 

/>oo 

Xi = c / n 1 (E 1 ,T)<r ( i 1 ' a) {E 1 )dE 1 

" ^WfWTT) af ( ± ^) f exp (-§ " 2 "') S ' lE " ,dE ° (6 - 2) 

with the thermal photon density 



1 1 E 2 

n 7 (E 7 , T) dEl = - 2W ^ e ^ {E J kT) _ 1 dEL, ( 1 
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(see e.g. |33), the excitation energy E x/ i and the spin 7j of the z-th state, the Sommerfeld parameter 77, and the 
astrophysical S-factor Si for the transition to the z-th state in 20 Ne. The derivation of Eq. I|6.2|) uses only the time 
reversal symmetry for the cross sections of transitions between individual states 0, 0] . Note that the energies are 
are related by 

E 7 = Q + E a - E Xti (6.4) 

with E a in the center-of-mass system. This means that the required photon energy E 7 for the photodisintegration of 
20 Ne is reduced by the excitation energy E x> i of the i-th state in 20 Ne. Consequently, the reaction rate A; of the i-th 

state is enhanced by a factor exp (^+-}$f-^j because of the enhanced photon density at lower photon energies. 

After summing all the rates according Eq. (|6.2|> with their proper thermal weights proportional to the Boltzmann 

factor cxp ^— ^f-^j under stellar conditions one obtains a stellar reaction rate 

x *~i^jKcf ^ exp {-w) f exp (-§ - H S ^ E " ] dE ° (6 - 5) 



with the total astrophysical S-factor S t ot = Si &i for the 16 O(a,7) 20 Ne capture reaction. The comparison of Eq. I|6.5|) 
with the usual definition of the reaction rate of the capture reaction leads directly to the detailed balance result in 
Eq. By the way, the analysis of the integrand in Eq. (|6.5|l together with Eq. (|6.4|) shows that the most effective 

energy for the inverse photodisintegration reaction is located at E x ( 20 Ne) = Q + Eq where Eq is the most effective 
energy of the capture reaction as defined in Sect. [I] Further details on the Gamow window for (7, a) reactions see [35j. 

Although their thermal occupation probabilities are extremely low because of the Boltzmann factor, excited states 
play a significant role in the photodisintegration process because the Boltzmann factor is exactly balanced by an 
enhancement factor because of the lower required photon energy according to Eqs. I|6.3I) and l|6.4l) . In the case of 
the 20 Ne(7,a) 16 O photodisintegration reaction, the first excited state in 20 Ne contributes with 78% to the photo- 
disintegration rate A* because of the 78% branching to this state in the 16 O(a,7) 20 Ne capture reaction although 
the thermal occupation probability e.g. at Tg = 0.2 is below 10 -40 ! But there is no significant contribution of the 
first excited state in 16 neither to the 16 O(a,7) 20 Ne capture reaction nor to the 20 Ne(7,a) 16 O photodisintegration 
reaction because there is no compensation in the cross section or reaction rate for the suppression factor from the 
Boltzmann statistics. 

Experimentally, photodisintegration rates in the laboratory can be measured using quasi-thermal photon spectra 
[3311361 ]. Contrary to the capture rate (see Sect. IIVB|l . the photodisintegration rates A* under stellar conditions and 
Ai a b in the laboratory are not identical: A* ^ Ai a b ! The laboratory experiment provides Ai a b with the target 20 Ne in 
its ground state which is given by Eq. I|6.2[l for i — g.s.; compared to Eq. I|6.5|l one finds a ratio of 

A*(T) ~ S tot (Eo) 1 ' 

around T 9 = 0.2 for the 20 Ne(7,a) 16 O reaction. E is the most effective energy of the capture reaction at temperature 
T, see Sect.Q] Provided that the branching ratio Sg.s./Stot is experimentally known, the stellar rate A* can easily be 
derived from the experimental result Ai a b- Otherwise, theoretical models have to be used to determine A* from Ai a b- 
In any case, experimental data may provide helpful information to restrict theoretical calculations of the laboratory 
cross section or rate [37T |. Further information on photodisintegration experiments for the 20 Ne(7,a) 16 O reaction is 
given in [32l| . 



VII. SUMMARY AND CONCLUSIONS 



The cross section of the 16 O(a,7) 20 Ne capture reaction has been calculated using a simple two-body DC model 
together with systematic folding potentials. The adjustment of the potential strength and the reduced transition 
probabilities to the relevant experimental quantitites enables the calculation of all relevant El, E2, and Ml transitions 
between incoming partial waves with < L < 6 and all bound states (J 71 " — + , E x — 0; 2 + , 1634keV; and 4 + , 
4248 keV) in 20 Ne without further free parameters. 

The astrophysical S-factor at energies below 1 MeV is somewhat higher than in previous calculations 0> 0, IS S j 
whereas the energy dependence and the branching ratio to the final states is in good agreement. Consequently, the 
astrophysical reaction rate is enhanced compared to a recent compilation Q] at temperatures below Tg — 0.2 where 
the direct capture contribution is dominating. 

The astrophysical reaction rate of the inverse 20 Ne(7,a) 16 O photodisintegration reaction has a significant contri- 
bution of the thermally excited first excited state in 20 Ne although the thermal occupation probability seems to be 
negligible. The widely used detailed balance theorem remains valid also for the case of high-lying first excited states. 
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